/******************************************************************************/
/*           GOLEM - Multiphysics of faulted geothermal reservoirs            */
/*                                                                            */
/*          Copyright (C) 2017 by Antoine B. Jacquey and Mauro Cacace         */
/*             GFZ Potsdam, German Research Centre for Geosciences            */
/*                                                                            */
/*    This program is free software: you can redistribute it and/or modify    */
/*    it under the terms of the GNU General Public License as published by    */
/*      the Free Software Foundation, either version 3 of the License, or     */
/*                     (at your option) any later version.                    */
/*                                                                            */
/*       This program is distributed in the hope that it will be useful,      */
/*       but WITHOUT ANY WARRANTY; without even the implied warranty of       */
/*        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the       */
/*                GNU General Public License for more details.                */
/*                                                                            */
/*      You should have received a copy of the GNU General Public License     */
/*    along with this program.  If not, see <http://www.gnu.org/licenses/>    */
/******************************************************************************/

#include "GolemHeatFlowBC.h"
#include "Function.h"
#include "GolemScaling.h"

registerMooseObject("GolemApp", GolemHeatFlowBC);

template <>
InputParameters
validParams<GolemHeatFlowBC>()
{
  InputParameters params = validParams<NeumannBC>();
  params.addParam<FunctionName>("function", "The function of the heat flow value.");
  params.addParam<UserObjectName>("scaling_uo", "The name of the scaling user object.");
  return params;
}

GolemHeatFlowBC::GolemHeatFlowBC(const InputParameters & parameters)
  : NeumannBC(parameters),
    _has_scaled_properties(isParamValid("scaling_uo") ? true : false),
    _function(isParamValid("function") ? &getFunction("function") : NULL),
    _scaling_uo(_has_scaled_properties ? &getUserObject<GolemScaling>("scaling_uo") : NULL)
{
}

Real
GolemHeatFlowBC::computeQpResidual()
{
  if (_has_scaled_properties)
  {
    if (isParamValid("function"))
      _scaled_value = _function->value(_t, Point()) / _scaling_uo->_s_heat_flow;
    else
      _scaled_value = _value / _scaling_uo->_s_heat_flow;
  }
  else
  {
    if (isParamValid("function"))
      _scaled_value = _function->value(_t, Point());
    else
      _scaled_value = _value;
  }

  return -_test[_i][_qp] * _scaled_value;
}